The absolute frequency of the 87 Sr optical clock transition 
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The absolute frequency of the 1 S'o— 3 i-*o clock transition of 87 Sr has been measured to be 
429 228 004 229 873.65 (37) Hz using lattice-confined atoms, where the fractional uncer- 
tainty of 8.6x 10~ 16 represents one of the most accurate measurements of an atomic transi- 
tion frequency to date. After a detailed study of systematic effects, which reduced the total 
systematic uncertainty of the Sr lattice clock to 1.5 x 10~ 16 , the clock frequency is measured 
against a hydrogen maser which is simultaneously calibrated to the US primary frequency 
standard, the NIST Cs fountain clock, NIST-F1. The comparison is made possible using a 
femtosecond laser based optical frequency comb to phase coherently connect the optical and 
microwave spectral regions and by a 3.5 km fiber transfer scheme to compare the remotely 
located clock signals. 



I. INTRODUCTION 



In recent years optical atomic clocks have made great strides, with dramatic improvements 
demonstrated in both stability and accuracy, and have now surpassed the performance of the best 
microwave standards HI 13- Optical clock candidates are being investigated by a variety of groups 
using a number of different atomic transitions in trapped ions [|2l[3l|4l[5l[6||, trapped neutral atoms, 
and freely expanding neutral atoms JB HI 13 OH El [HI H3- As the best optical standards 
now support an accuracy surpassing that of the Cs primary standards (3.3 x 10~ 16 ) lfT4l [T51 [T6l , 
it becomes imperative to directly compare these optical standards against each other HI 13 to 
evaluate them at the lowest possible level of uncertainties. Nevertheless, it still remains important 
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that these optical standards are evaluated by the mature primary Cs standards for multiple reasons. 

First, the accuracy of frequency standards is ultimately defined by the Cs clock under the cur- 
rent realization of the Si-second. Additionally, over the years a remarkable infrastructure has been 
developed to support the transfer of Cs standards for international intercomparisons, and the pri- 
mary frequency standards at multiple national labs all agree within their stated uncertainties ifTTI . 
While fiber networks [fT8l [T9l l20l ETI now provide the most precise frequency distribution links 
between optical clocks located near each other (for example within 100 km), for intercontinental 
comparisons optical clocks need to be measured relative to Cs standards. In fact, recent intercom- 
parisons of Sr clocks among three laboratories at JILA, SYRTE, and University of Tokyo [22J 
have reached an agreement at 1 x 10~ 15 , approaching the Cs limit. This has firmly established the 
Sr lattice clock standard as the best agreed-upon optical clock frequency to date, and second only 
to the Cs standard itself among all atomic clocks. 

Second, an important application of highly accurate atomic clocks is the test of fundamen- 
tal laws of nature with high precision. For example, atomic clocks are placing increasingly 
tighter constraints on possible time-dependent variations of fundamental constants such as the 
fine-structure constant (a) and the electron-proton mass ratio (ji) ll2ll2^l231 IMl231 IM l27ll28ll29ll . 
These measurements are made by comparing atomic transition frequencies among a diverse set of 
atomic species, helping reduce systematic effects. For example, an optical clock transition fre- 
quency is generally sensitive to variations of a, with different atoms having different sensitivities 
1501 . Sr in fact has a rather low sensitivity. The Cs standard on the other hand is based on a hyper- 
fine transition and is sensitive to variations in both a and /i. Thus measurement of the frequency 
ratio of Sr and Cs over the course of a year limits not only the possible linear drift of these con- 
stants but also constrains possible coupling between fundamental constants and the gravitational 
potential, which would signal a violation of local position invariance [|22ll29l[3Tl[32l . 

In recent years, the most accurate absolute frequency measurements were performed using 
single trapped ions. These systems benefit from the insensitivity of the ions to external pertur- 
bations, and using Hg + ions a frequency uncertainty of 9.1 xlO -16 has been achieved. Large 
ensembles of neutral atoms offer high measurement signal to noise ratios, however, neutral atom 
systems have typically been limited by motional effects. By confining the atoms in an optical 
lattice BH [H |H these effects are greatly reduced, as the atoms can be trapped in the Lamb-Dicke 
regime, where both Doppler and photon-recoil related effects are suppressed. One such system is 
the 87 Sr (5s 2 yS — (5s5p) 3 P transition, which is currently being pursued by a number of groups 
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worldwide EES 133 EH- 

In this paper we report on the absolute frequency measurement of the 87 Sr 1 S — 3 P clock 
transition. The absolute frequency is measured using a femtosecond laser based octave-spanning 
optical frequency comb to compare the 87 Sr optical transition frequency to a hydrogen maser, 
which is simultaneously calibrated to the NIST fountain primary frequency standard, NIST-F1. 
To remotely link the Sr standard, which is located at JILA on the University of Colorado campus, 
to the NIST-F1 Cs clock, located at the NIST Boulder laboratories, a 3.5 km optical fiber link is 
used to transfer the H-maser reference signal [|36l [37). In addition to demonstrating one of the 
most accurate measurements of an optical transition frequency to date, the agreement of this result 
with previous measurements both at JILA and around the world demonstrates the robustness and 
reproducibility of strontium as a frequency standard, and as a future candidate for the possible 
redefinition of the SI second. 



II. EXPERIMENTAL SETUP 



The frequency standard uses lattice-confined 87 Sr atoms with nuclear spin / = 9/2. Although 
the 87 Sr apparatus has been previously described elsewhere [[8l [33), here we summarize the ex- 
perimental details most relevant to this work. To measure the frequency of the clock transition, 
87 Sr atoms are first trapped and cooled to mK temperatures in a magneto-optical trap (MOT) op- 
erated on the 1 Sq— 1 P\ strong 461 nm cycling transition (see Fig. [T^ for a diagram of relevant 
energy levels). The atoms are then transferred to a second stage 689 nm MOT for further cooling. 
This dual-frequency MOT uses narrow line cooling 11381139) , resulting in final temperatures of ~1 
fjK. During the cooling process, a one-dimensional optical lattice is superimposed in the nearly 
vertical direction. After the second MOT stage, the MOT optical beams and the inhomogeneous 
magnetic field are turned off, leaving ~10 4 atoms at 2.5 /xK trapped in the optical lattice. The 
optical lattice is created using a retro-reflected laser beam and is operated near a laser frequency 
where the polarizability of the 1 Sq and 3 P states are identical for the lattice field II401I4D . For this 
work, the lattice is operated at a trap depth of Ut = 35E rec , where E rec = h 2 k 2 /2m is the lattice 
photon recoil energy and k = 2n/\is the wavevector of the lattice light. At this lattice depth the 
atoms are longitudinally confined in the Lamb-Dicke regime and in the resolved sideband limit 
j|42~). Spectroscopy is performed by aligning the probe laser precisely along the axis of the lattice 
standing wave, and the atoms are probed free of recoil or motional effects. The vertical orienta- 
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tion of the lattice breaks the energy degeneracy between lattice sites, strongly prohibiting atomic 
tunneling ll43l 

Before performing spectroscopy, the atoms are first optically pumped to the stretched \F = 
9/2, nip = ±9/2) states with the use of a weak optical beam resonant with the l So(F = 9/2)- 3 Pi 
(F = 111) transition. Here F = I + J is the total angular momentum, with / the nuclear spin and 
J the total electron angular momentum. The beam used for optical pumping is aligned collinear 
with the lattice, and is linearly polarized along the lattice polarization axis. The optical pumping 
is performed with a small magnetic bias field (~3 /iT), which is also oriented along the lattice 
polarization. After optical pumping, spectroscopy is performed on the 1 S - 3 P clock transition 
from the two spin sublevels. The clock transition, which has a theoretical natural linewidth of 
~1 mHz [|4^1l44ll45ll46ll . is interrogated using a diode laser at 698 nm, which is prestabilized by 
locking it to a high-finesse ultrastable cavity, resulting in a laser optical linewidth below 1 Hz [|47l . 
The probe beam is coaligned and copolarized with the optical lattice. To ensure that the stretched 
states are well resolved, the spectroscopy is performed under a magnetic bias field of 25 /iT, which 
results in a ~250 Hz separation between the two 7r-transitions excited during the spectroscopy. 

Spectroscopy is performed using an 80-ms Rabi pulse, which when on resonance transfers a 
fraction of the atoms into the 3 P state. After applying the clock pulse, atoms remaining in the 1 So 
ground state are detected by measuring fluorescence on the strong 1 So— l P\ transition. The length 
of the pulse is long enough to measure both the population in the 1 S state as well as to heat these 
atoms out of the trap. The population in the 3 P state is then measured by first pumping the atoms 
back to the 1 S state through the intermediate (5s5p) 3 P —(5sQs) 3 Si and (5s5p) 3 P 2 — (5s6s) 3 Si 
states and then by again measuring the fluorescence on the 1 So— l P\ transition. Combining these 
two measurements gives a normalized excitation fraction insensitive to atomic number fluctuations 
from shot to shot. A typical spectrum is shown in Fig. [1]). The Fourier-limited linewidth of the 
transition is 10 Hz, much less than the 250 Hz separation between the peaks, which makes the lines 
well resolved and also reduces potential line pulling effects due to any residual population left in 
other spin states by imperfect optical pumping. We note that while our optical local oscillator 
supports recovery of < 2 Hz spectroscopic linewidths [|48l . we find it more robust to run the clock 
transition with a 10 Hz Fourier- limited spectral linewidth. 

To stabilize the optical local oscillator used for spectroscopy to the atomic transition, we use 
both stretched states. Using two time-multiplexed independent servos, we lock the laser frequency 
to the center of each transition. This is done by sampling the full width half maximum (FWHM) 
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of each transition (labeled f\ st i oc k in Fig. [lj)). The average of the two line centers gives the 
center frequency of the clock transition. The cavity- stabilized local oscillator, in combination 
with this frequency f center, is in turn used to phase-lock a self-referenced octave- spanning optical 
frequency comb. The Sr-referenced repetition frequency of the comb is then counted with a H- 
maser located at NIST. A schematic of this locking setup is shown in Fig. [Tb. Determination of 
the center frequency requires four experimental cycles, two for each of the mp = ±9/2 transitions 
since a new atomic sample is reloaded for each lock point. The length of each experimental cycle 
is ps 1.1s. After first probing the n transition for the m F = —9/2 transition, we then probe the 
transition for the m F = +9/2 state. The digital servo operates via standard modulation techniques. 
A linear feedback is also implemented to compensate for the drift of the high-finesse cavity used 
to prestabilize the clock laser. A second integration stage in the laser-atom feedback loop is used 
to calculate this feedback value (labeled /2nd lock) in Fig. [lb. As shown in Fig. fTti, using this 
approach limits the residual drifts compensated for by the first servo integrator to typically < 1 
mHz/s. 



III. SYSTEMATIC SHIFTS OF THE STRONTIUM CLOCK 

We have recently evaluated the systematic shifts of the strontium clock at the 1 x 1CT 16 level 
HI, and in Table [I] the important systematic shifts to the absolute frequency are shown. Although 
a detailed description can be found in HI, here we summarize these shifts. The evaluation of the 
systematic uncertainty is performed using the remotely located calcium optical standard at NIST 
fl49), which is linked to JILA via a phase coherent optical fiber link [18]. The Sr-Ca comparison 
has a 1 s stability of 2x 1CT 15 which averages down to below 3 x 1CT 16 after 200 s. To measure the 
Sr systematics an interleaved scheme is used where the Sr parameter of interest is varied between 
two different settings every 100 s, while the Ca standard remains locked. Pairs of such data are 
then used to determine the frequency shift, and many pairs are averaged in order to reach below 
the 1(T 16 level. 

As shown in Table [j] besides the correction that arises from the maser calibration, the dominant 
shift for the Sr clock is the black-body radiation (BBR)-induced shift. To determine this shift, 
the temperature of the Sr vacuum chamber is continuously monitored during the course of the 
absolute frequency measurement at four separate locations. During the measurement the average 
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TABLE I: Frequency corrections and their associated uncertainties for the clock transition in units of 10 
fractional frequency, and with vq =429 228 004 229 800 Hz. The maser correction uncertainty includes 
both Sr/H-Maser comparison, as well as the Cs clock uncertainty. 



temperature of the chamber is 295(1) K, and the corresponding BBR effect: 



T ^ 



Svbbr = -2.354(32) J Hz (1) 

gives a frequency shift of 54(l)xlCT 16 . Higher-order multipoles are suppressed by a 2 and are 
negligible at this level. The given uncertainty in the BBR shift includes the error due to the 
chamber temperature, as well as the theoretical uncertainty in the polarizability ll50l . 

For the duration of the experiment, the lattice laser is phase locked to the same optical frequency 
comb used to count the Sr beat, and the wavelength is simultaneously monitored on a wavemeter 
to ensure it does not mode-hop. The lattice is operated at a frequency of 368554.36(3) GHz 03 [EL 
slightly away from the state-insensitive lattice frequency. Including nuclear spin effects, the light 
shift due to the linearly polarized lattice can be expressed as [|44l : 

5u s « -(Ak s - Ak t F(F + 1))-^- - (A/m^ cos(^) + A« T 3m|)—^- (2) 

^rec -^rec 

where Ak s ' t ' v is the frequency shift coefficient due to the differential polarizability (scalar, tensor, 
and vector) between the ground and excited clock states, £ is the degree of ellipticity of the beam 
(with £ = for 7r -polarized light), and ip is the angle between the lattice propagation direction 
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and the bias magnetic field (~ 7r/2 for our setup). For the linearly polarized lattice configura- 
tion in our setup, the vector light shift is minimized; furthermore since the Sr clock is operated 
using both the m F = ± 9/2 states, the antisymmetric m F dependence averages away this vector 
shift. The effect of the tensor light shift for a given \t71f\ state introduces a polarization dependent 
offset to the state-insensitive lattice frequency. Experimentally, this Stark cancelation frequency 
for the m F = ±9/2 state has been determined to be 368554.68(17) GHz [1J leading to a shift of 
-6.5(5) x 10 -16 to the absolute clock frequency for the lattice depth and frequency used during this 
measurement. For our operating conditions, hyperpolarizability effects are more than an order of 
magnitude smaller [|9l|4Tj. The ground and excited clock states have different polarizabilities at the 
clock transition frequency, and imperfect alignment between the clock laser and the lattice beam 
can lead to inhomogeneous Rabi frequencies in the transverse direction requiring an increase in 
the clock transition probe power. However, given the small saturation intensity of the clock probe 
beam, the ac Stark shift introduced by the clock laser during spectroscopy is small, and has been 
experimentally measured to be 0. 15(10) xlCT 16 . Stark shifts from laser beams not used during 
spectroscopy, for example those used for cooling and trapping during the MOT phase and for fluo- 
rescence detection after spectroscopy, are eliminated through the use of acousto-optic modulators 
(AOM) in series with mechanical shutters which block these beams during spectroscopy. In addi- 
tion, the vacuum chamber is covered with an opaque cloth to prevent any stray light from entering 
the chamber. 

For each experimental cycle sequence, the total atom number is recorded, allowing a point-for- 
point correction of the density shift. The value given in Table [I] is the average density correction. 
At the FWHM of the spectroscopic signal where the probe laser is locked, the excitation frac- 
tion in each stretched state is 15(2)%. This excitation fraction and our operating density of ~ 
4x 10 10 cm~ 3 gives a frequency correction of 3.8(5)x 1CT 16 flUED. 

The Zeeman shift of the transition frequency is given in Hz by [|44l 

bv B « -5gm F fi B + (3B 2 « -1.084(4) x 10 6 m F B - 5.8(8) x 10~ 8 B 2 (3) 

where /i = ^b/K with fx B the Bohr magneton, and 5g the differential Lande g-factor between 
the ground and excited clock states. The l st -order Zeeman shift is experimentally measured in 
P4ll . and the 2™ d -order Zeeman shift is experimentally measured in 0], consistent with other 
measurements ll52~ll . By measuring the average frequency of the stretched states at a small bias 
field, the 1 st order Zeeman effect is averaged away due to the opposite linear dependence of the 
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shift on the nip = ±9/2 states, and the experimentally measured value for the shift is consistent 
with zero. The bias field of 25 fiT used during spectroscopy is large enough such that the spin 
states are well resolved, reducing line pulling effects due to residual populations in other states, 
yet small enough such that the 2 nd order Zeeman shift is negligible, with a value of 0.36(4) x 10~ 17 
for our bias field. 

By operating in the Lamb-Dicke regime, l si order Doppler shifts are minimized. However, 
driven motion can also cause frequency shifts due to shaking of the lattice beams, or due to relative 
motion between the lattice and the probe beam. To minimize vibrations, the optics table is floated 
using standard pneumatic compressed air legs, and we estimate the effect of I s * order Doppler 
shifts to be below 1CT 18 . Switching magnetic fields can also induce vibrations, however, our 
quadrupole trap is switched off more than 100 ms before spectroscopy. Furthermore, 2 nd -order 
Doppler effects from residual thermal motion are negligible (< 1CT 18 ), given the T = 2.5 fiK 
temperature of the trapped Sr atoms. 

The digital servos used to steer the spectroscopy laser to the atomic transition are another 
potential source of frequency offsets. The dominant cause of servo error is insufficient feedback 
gain to compensate for the linear drift of the high-finesse reference cavity. The second integration 
step as described in the experimental setup section reduces this effect. By analyzing our servo 
record we conservatively estimate this effect to be < 5 x 10~ 17 . In conclusion, with the exception 
of the BBR-induced shift, all of the systematics discussed above are limited only by statistical 
uncertainty. 

IV. FIBER TRANSFER BETWEEN JILA AND NIST 

The strontium experiment is located at JILA, on the University of Colorado, Boulder campus 
and is linked to NIST Boulder Laboratories by a 3.5 km optical fiber network. To measure the 
absolute frequency of the transition, the tenth harmonic of the repetition rate of the frequency 
comb (which is phase-locked to the Sr clock laser and located at JILA), is beat against a ~950 
MHz signal originating from NIST. A schematic of the transfer scheme from NIST is shown 
in Fig. [2} A 5 MHz signal from a hydrogen maser is distributed along a ~300 m cable to a 
distribution room where it is first frequency doubled and then used as a reference to stabilize an 
RF synthesizer operating near 950 MHz. In Fig. [2}), the stability of the H-maser as measured by 
Cs is shown, where the total length of the Sr absolute frequency measurement is indicated with a 
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dotted line. The RF signal generated by the synthesizer is then used to modulate the amplitude of 
a 1320 nm diode laser which is then transferred to JILA via the 3.5 km optical fiber link between 
the two labs 11361 |37ll . The microwave phase of the fiber link is actively stabilized using a fiber 
stretcher to control the group delay between NIST and JILA ll37ll . The limited dynamic range 
of the fiber stretcher necessitates a periodic change of the transfer frequency to relock the fiber 
transfer system for our 3.5 km link. Typically the dynamic range is sufficient to stabilize drifts for 
roughly 30-60 minute intervals, after which it must be unlocked and reset, leading to a dead time 
in the measurement of ~1 minute. 

The transfer of the microwave signal between the NIST H-maser and JILA can potentially in- 
troduce a number of systematic frequency shifts and uncertainties. The majority of these arise 
from temperature-driven fluctuations during the course of the measurement. The microwave sig- 
nal is transferred between the maser and the RF synthesizer using 300 m of cable after which it 
goes through a series of distribution amplifiers and a frequency doubler. All of the microwave 
electronics, as well as the cable used to transfer the signal, are sensitive to temperature-driven 
phase excursions. In order to correct for these effects, the temperature in the room is monitored 
continuously during the course of the experiment. In addition, the RF synthesizer is placed in 
a temperature- stabilized, thermally insulated box, and the temperature in the box is also contin- 
uously monitored. The temperature coefficient of the synthesizer is independently measured by 
applying a temperature ramp to the box while counting the frequency of the synthesizer relative to 
a second frequency stable synthesizer. The synthesizer is found to have a temperature coefficient 
of -3.6 ps/K, corresponding to a fractional frequency change of -1 x 1CT 15 for a temperature ramp 
of 1 K/h. This temperature coefficient is used to make rolling frequency corrections during the 
absolute frequency measurement. 

To test the performance of the microwave electronics used for the modulation of the transfer 
laser as well as the fiber noise cancelation, an out-of-loop measurement is performed by detecting 
the heterodyne beat between the resulting transfer signal and the RF synthesizer. In Fig. [2£ the 
Allan deviation is shown for this measurement, demonstrating that the fractional frequency insta- 
bility due to these transfer components is 1 x 1CT 14 from 1-10 s, and averages down to < 1CT 17 at 
10 5 s. By correlating the out-of-loop measurement with temperature fluctuations in the distribu- 
tion room during the course of the measurement, a temperature coefficient of 4.4 x lCT 16 /(K7h) is 
found for the microwave electronics. 

The distribution amplifiers, frequency doubler, and cable used to transfer the signal within the 
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distribution room are also tested by comparing the 10 MHz signal used to stabilize the synthesizer 
with a signal split off before the distribution amplifiers. This measurement determined a coefficient 
of ~ 7 x 10~ 16 /(K/h). 

The absolute frequency measurement, as discussed below in Section V, is recorded over 50 
continuous hours. During the course of the absolute frequency measurement the insulated box used 
to house the RF synthesizer maintained an average temperature of 292.2(2) K, with a maximum 
slope of < 0.1 K/h. In Fig.[2]i the black trace shows the resulting fractional frequency correction to 
the Sr frequency due to temperature-driven frequency fluctuations of the RF synthesizer. The total 
fractional correction to the Sr frequency due to the microwave electronics during the measurement, 
as well as the distribution amplifiers and cables, is shown by the grey trace in Fig. |2ji. During 
hour 37 (133 xlO 3 s) of the measurement, a temperature ramp began in the distribution room, 
leading to a large slope in the temperature during hour 37 and during hour 50 (180xl0 3 s) as 
the temperature restabilized. However, this transient affected only a small fraction of the data. 
Using the measured temperature coefficients, a rolling correction is made to all of the measured 
frequencies. The average frequency correction during the course of the measurement is 9 x 10~ 17 , 
with an uncertainty of 1 x 10~ 17 . These corrections do not influence the statistics of the final 
absolute frequency measurement. In other words, the final mean frequency and standard error are 
the same with and without these corrections. 

V. FREQUENCY MEASUREMENT RESULTS 

As shown in Fig. [2J5, the Allan deviation of the H-maser averages down as ~ 3 x 10~ 13 / yfr. To 
measure the Sr absolute frequency to below 10~ 15 , the measurement is performed for 50 continu- 
ous hours. The largest frequency correction to the measured 87 Sr frequency is the calibration offset 
of the H-maser frequency. The H-maser is simultaneously counted against the Cs standard during 
the duration of the measurement and the resulting frequency correction to the 87 Sr/H-maser com- 
parison is -439.37(85)x 10~ 15 , where the uncertainty includes both an uncertainty of 0.6xl0~ 15 
due to the Cs standard, as well as an uncertainty of 0.6 x 10~ 15 from dead time in the Sr/H-maser 
measurement. An additional frequency correction is the gravitational shift due to the difference 
in elevation between the Cs laboratory at NIST and the Sr laboratory at JILA. The difference in 
elevation between the two labs, which has been determined using GPS receivers located in each 
building to be 1 1.3(2) m, gives a frequency shift of 12.5(1.0)x 10~ 16 , where the given uncertainty 
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includes the uncertainty in the elevation as well as the uncertainty due to the geoid correction ll53l . 

In Fig. [3j the 50 hour counting record is shown for the Sr frequency, with a 30 s gate time. The 
frequency shown includes only the correction due to the maser, and is plotted with an offset fre- 
quency of uq = 429 228 004 229 800 Hz. The frequency excursions and gaps seen in Fig. [3^ occur 
when the Sr system is unlocked, which happens when either the frequency comb comes unlocked, 
or when the probe laser is not locked to the atomic signal. During the course of the measurement 
the lattice intensity and frequency, all laser locks, and the temperature at both JILA and in the dis- 
tribution room in NIST are continuously monitored and recorded. In Fig. [3}), data corresponding to 
times when any lasers are unlocked, including times when the spectroscopy laser is not locked to 
the atoms, and times when the lattice laser intensity or frequency is incorrect, have been removed. 
In Fig. [4^, a histogram of this final counting record is shown, demonstrating the gaussian statistics 
of the measurement. The mean value (relative to vq) of the measured frequency is 70.88(35) Hz. 
In Fig. |4j5 the Total deviation of the frequency measurement is shown. The Total deviation [|54l is 
similar to the Allan deviation, however it is better at predicting the long-term fractional frequency 
instability. The 1-s stability of the H-maser used for the measurement is 1.5xl0~ 13 . However, 
from a fit to the Total deviation, we find a 1-s stability of o\ a = 2.64(8) x 10 -13 for the counting 
record, which is limited by counter noise, and averages down as a\ s /r 0-48 ^. Extrapolating to the 
full length of the data set (excluding dead time in the measurement) gives a statistical uncertainty 
of 8xlCT 16 . The frequency uncertainty of Sr is low enough such that this uncertainty is domi- 
nated by the performance of the maser (see Fig. |4j)), which is included in the maser calibration 
uncertainty given in Table [IJ 

Including the uncertainty of the H-maser as well as the strontium systematics described in Sec- 
tion III gives a final frequency of 429 228 004 229 873.65 (37) Hz, where the dominant uncertainty 
is due to the H-maser calibration. In Fig. [5j this measurement is compared with previous Sr fre- 
quency measurements by this group J8l |33l [55l |56]|, as well as by the Paris [33J and Tokyo 
||34l groups. As shown in the figure, the agreement between international measurements of the 
Sr frequency is excellent, with the most recent measurements in agreement below the 10~ 15 level, 
making the Sr clock transition the best agreed upon optical frequency standard to date. The high 
level of agreement enabled a recent analysis of this combined data set that constrains the coupling 
of fundamental constants to the gravitational potential as well as their drifts [|22l . 
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VI. CONCLUSION 

In conclusion, we have made an accurate measurement of the 1 S — 3 P clock transition in 
fermionic strontium, where the final fractional uncertainty of 8.6 x 10~ 16 is limited primarily by the 
performance of the intermediate hydrogen maser used to compare the Sr standard to the NIST-F1 
Cs fountain clock. This experiment represents one of the most accurate measurements of an optical 
frequency to date, and the excellent agreement with previous measurements makes strontium an 
excellent candidate for a possible redefinition of the SI second in the future. In addition, the 
combined frequency measurements of 87 Sr performed worldwide, as well as future measurements 
of frequency ratios with other optical standards can be used to search for time-dependent frequency 
changes which constrain variations of fundamental constants ||22]|. 
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FIG. 1: Experimental Setup, (a) Relevant energy levels for 87 Sr used for the optical lattice clock. Transi- 
tions at 461 nm and 689 nm are used in two-stage cooling and trapping of the Sr atoms. The clock transition 
is at 698 nm. Lasers at 679 and 707 nm provide necessary repumping from metastable states, (b) To operate 
the clock, ultracold 87 Sr atoms are first optically pumped to the \F = 9/2, m/ = ±9/2) states. The clock 
center frequency (f center) is found by locking the probe laser frequency to both peaks successively and tak- 
ing their average. The laser is locked to the center of each transition by sampling their FWHM as illustrated 
in the figure by dots (fist lock)- (c) Schematic of the setup used for locking the optical local oscillator to the 
87 Sr transition. The clock transition is probed using a diode laser (DL) at A = 698 nm which is locked to an 
ultrastable, high finesse optical cavity. The laser beam is used to interrogate the Sr atoms and is transferred 
to the atoms using an optical fiber with active fiber noise cancelation. To steer the frequency of the laser for 
the lock to the Sr resonance, an acousto-optic modulator (AOM) is used to introduce a frequency offset be- 
tween the cavity and the atoms. The frequency offset is steered to the lock points (fist lock)- The frequency 
offset also includes a linear feedback value (/2nd lock) to compensate for the linear drift of the high finesse 
cavity. The cavity-stabilized clock laser is also used to phase-lock a self -referenced octave-spanning optical 
frequency comb, in combination with atomic resonance information contained in f cen ter and f^nd lock- The 
Sr-referenced repetition frequency of the comb (/ rep ) is then counted relative to a H-maser located at NIST 
(fmaser)- (d) Sample data showing the in-loop atom lock for 5000 s of data taken during the measurement 
of the absolute frequency. The fit gives a residual linear fractional frequency drift of < 2xl0 _19 /s. 
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FIG. 2: Clock signal transfer between NIST and JILA. (a) Schematic of the setup used to transfer the 
hydrogen maser signal from NIST to JILA. A 5 MHz signal from the H-maser, which is simultaneously 
counted against the NIST-F1 Cs standard, is distributed through a ~300 m cable to a distribution amplifier 
(DA). After the distribution amplifier it is actively frequency doubled (2f) and sent through a second dis- 
tribution amplifier. The resulting 10 MHz signal is used to reference an RF Synthesizer operating at ~950 
MHz. The synthesizer in turn modulates the amplitude (AM) of a 1320 nm laser (DL) which is transferred 
to JILA through a 3.5 km fiber link. Noise from the fiber link is canceled with a fiber stretcher to actively 
stabilize the microwave phase using a retroreflection of the beam sent back from JILA [37]. (b) Typical 
Allan deviation of the H-maser used for the Sr absolute frequency measurement, with the total duration of 
the measurement represented by a dotted line, (c) Out of loop measurement of the stability of the microwave 
electronics used for transfer and fiber noise cancelation. The fit to the line gives a 1-s Allan deviation of 
1.08(1) xlO -13 with a slope of -0.889(4). The bump at 10000 s is indicative of temperature fluctuations in 
the distribution room during the out of loop measurement. The bump at 10 s is due to low pass filtering of 
the phase measurement, (d) Frequency correction due to temperature fluctuations in the distribution room 
(grey curve) and fluctuations in the temperature-stabilized box used to house the RF synthesizer (black 
curve) during the course of the measurement. 
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FIG. 3: Absolute frequency measurements of the 1 So— 3 Po clock transition, (a) Counting record showing 
all of the data taken over a 50 hour period. Each point corresponds to a 30 s average, and the overall offset 
is vq = 429 228 004 229 800 Hz. (b) The counting record after removing points where the system is not 
locked. The mean value is 70.88(35) Hz. 
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FIG. 4: (a) Histogram of the frequency measurements shown in Fig. [3} including the maser correction. 
The dashed line is a gaussian fit to the data, the mean frequency is 70.88 Hz and is indicated by the black 
line, (b) Total Deviation of the frequency measurement for the Sr/H-maser comparison (circles), and the 
H-maser/Cs comparison (squares). The dotted line shows a fit of the Sr deviation to ar~ b , where a = 
2.64(8)xl0 -13 and b = 0.48(1), and the dotted line extends out to the full measurement time. For averaging 
times r > 10 4 s, the maser noise dominates both the Cs/H-maser and the Sr/H-maser measurement, and 
hence the maser uncertainty (6xl0 -16 , as described in the text) needs to be counted only once in the final 
Sr/Cs measurement uncertainty budget. 
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FIG. 5: Record of Sr absolute frequency measurements. Previous measurements by this group (circles) 
(U [33] SUES, as wen as tne Paris (triangle) 01351 and Tokyo (square) ll34ll groups are shown. The inset 
shows the high agreement of the most recent measurements which agree below the 1CP 15 level. The dashed 
line shows the weighted mean, v = 429 228 004 229 873.73 Hz of the combined data set. 



